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We review recent studies on exotic states at the Belle experiment. The results include: (1) The 
measurement of the cross sections of 77 — s- (0(j), 00, and (0(0 for masses that range from thresh- 
old to 4.0 GeV. In addition to signals from well established spin-zero and spin-two charmonium 
states, there are clear resonant structures below charmonium threshold, which have not been pre- 
viously observed. We report a spin-parity analysis for the new structures; (2) No Y(3872) signal 
is observed in rjJ/y/ or JXd mode in B decays. A narrow peak at 3823.5 MeV/c 2 (named as 
W2) to YXd with a significance of 4.2 standard deviations including systematic uncertainty is ob- 
served in ZP — ^ K^YXcll (3) The bottomonium states ht(\P), hb(2P) and T(1D) are observed in 
the reaction e + e~ —> k + k~ +X; (4) The observation of two narrow charged structures (named 
as Z 6 (10610) and Z 6 (10650)) in the mass spectra of the Tr^nS^n = 1,2,3) and n^h^mP) 
(m =1,2) pairs that are produced in association with a single charged pion in Y(5S) decays. 
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1. Introduction 

In hadronic physics, the best understood quark-antiquark systems are heavy quarkonia, i.e, 
cc or bb mesons. The discovery of the missing cc or bb states and the precise measurements of 
properties of the already observed ones are important. 

The QCD -motivated models predict the existence of hadrons of more complex structure than 
conventional mesons or baryons, such as hybrids, multiquark states of either molecular, tetraquark 
or hadrocharmonium configuration. As the conventional hadron spectrum is much cleaner than 
the dense spectrum of light states, exotic states containing cc or bb are expected to be identified 
more easily than the ones predicted in the light spectrum. Any resonance observed in addition to 
predicted multiplets might give a hint of such an exotic spectroscopy. 

Some of the recently observed charmonium-like or bottomonium-like XYZ states could be 
candidates for the exotic hadrons mentioned. However most of them still await confirmation or 
their properties need to be further studied before any decisive interpretation is made. Here, we 
review some recent results on the exotic hadrons from Belle experiment. 



2. Observation of new resonant structures in 77 — > Gxj), 00 and coco 

Recently a clear signal for a new state X(3915) — > (oJ/y ^ and evidence for another state 
X(4350) — > <j)J/y [Q] have been reported, thereby introducing new puzzles to charmonium or 
charmonium-like spectroscopy. It is natural to extend the above theoretical picture to similar states 
coupling to ft)0, (0(0 or 00. 

The measurements of the cross sections for 77—)- VV [||], where VV = 0)0,00 and (0(0, are 
based on an analysis of an 870 fb~ 1 data sample taken at or near the T(nS) (n = 1 , . . . , 5) resonances 
with the Belle detector operating at the KEKB asymmetric-energy e + e~ collider. 

After event selections, clear co or signal is observed. The magnitude of the vector sum of 
the final particles' transverse momenta in the e + e~ center-of-mass (CM.) frame, which 
approximates the transverse momentum of the two-photon-collision system, is used as a discrim- 
inating variable to separate signal from background. We obtain the number of VV events in each 
VV invariant mass bin by fitting the | £P*| distribution between zero and 0.9 GeV/c. The resulting 
VV invariant mass distributions are shown in Fig. [TJ where there are some obvious structures in the 
low VV invariant mass region. 
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Figure 1: The (a) 0)0, (b) and (c) (0(0 invariant mass distributions. The shaded histograms are from the 
corresponding normalized sidebands, which will be subtracted in calculating the final cross sections. 

Two-dimensional (2D) angular distributions are investigated to obtain the J p of the structures. 
In the process 77—^ VV, five angles are kinematically independent. We choose z, z* , z** , 0*, and 
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(j)** [Q] and use the transversity angle (fa) and polar-angle product (Tig) variables to analyze the 
angular distributions. They are defined as fa = \(j>* + <j)**\/27i, Yl e = [1 - (z*) 2 ][l - (z**) 2 ]. 

We obtain the number of signal events by fitting the distribution in each fa and Tig 

bin in the 2D space, which is divided into 4 x 4, 5 x 5, and 10 x 10 bins for co<p, (pip, and (0(0, 
respectively, for M(VV) < 2.8 GeV/c 2 , in some wider VV mass bins as shown in Fig. |[ The 
obtained 2D angular distribution data are fitted with the signal shapes from MC-simulated samples 
with different J p assumptions (0 + , 0~, 2 + , 2~). We find: (1) for coty: + (S-wave) or 2 + (S-wave) 
can describe data with % 2 /ndf = 1.1 or 1.2, while a mixture of + (S-wave) and 2+ (S-wave) 
describes data with % 2 /ndf = 0.9 (ndf is the number of degrees of freedom); (2) for <p<p : a mixture 
of + (S-wave) and 2~ (P-wave) describes data with % 2 /ndf = 1.3; and (3) for coco: a mixture of 
+ (S-wave) and 2 + (S-wave) describes data with % 2 /ndf = 1.3. 

The 77 — > VV cross sections are shown in Fig. The cross sections for different J p values as 
a function of M(VV) are also shown in Fig. |[ While there are substantial spin-zero components in 
all three modes, there are also significant spin-two components, at least in the <p<p and coco modes. 




3 

) (GeV/c 2 ) 



1 



(b) 



M{#) (GeV/c 2 ) 



data 

T 



Hi »i 



.5 2 2.6 




(<0 



5 2 2.5 3 3.5 4 
M((DG>) (GeV/c 2 ) 



Figure 2: The cross sections of yy —> (0(j) (a), 00 (b), and (0(0 (c) are shown as points with error bars. The 
cross sections for different J p values as a function of M(VV) are shown as the triangles and squares with 
error bars. The inset also shows the cross section on a semi-logarithmic scale. In the high energy region, 
the solid curve shows a fit to a Wyy n dependence for the cross section after the significant charmonium 
contributions (r/ c , XcO an d % c i) were excluded. 



The cross sections for yy — > co<p are much lower than the prediction of the q q 2 tetraquark 
model []|] of 1 nb, while the resonant structure in the 77 — » <p<p mode is nearly at the predicted 
position. However, the <pty cross section is an order of magnitude lower than the expectation in 
the tetraquark model. On the other hand, the t-channel factorization model [||] predicted that the 
<p<p cross sections vary between 0.001 nb and 0.05 nb in the mass region of 2.0 GeV/c 2 to 5.0 
GeV/c 2 , which are much lower than the experimental data. For 77 — > coco, the t-channel factoriza- 
tion model [||] predicted a broad structure between 1.8 GeV/c 2 and 3.0 GeV/c 2 with a peak cross 
section of 10-30 nb near 2.2 GeV/c 2 , while the one-pion-exchange model |7|] predicted an enhance- 
ment near threshold around 1.6 GeV/c 2 with a peak cross section of 13 nb using a preferred value 
of the slope parameter. Both the peak position and the peak height predicted in ^ and [Q] disagree 
with our measurements. 

3. Charmonium and charmonium-like states 

If X(3872) is a tetraquark state, then it has a C-odd parity (C = — ) partner, which can domi- 
nantly decay into J/yrj and % c \ y. B ± — > (J /yr\ (— > YY))^ an d — > (Xci (— > J/¥Y)y)K ± decay 
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modes are used in the search for C-odd partner of the X(3872) and other new narrow resonances. 
In all the decay modes, 7/y/ris reconstructed via e + e~ and . B candidates are identified using 
energy difference AE =Eg — E£ eam and beam-energy constrained mass M\, c = a/ (E? 7eam ) 2 — {p*b) 2 > 
where E£ eam is the beam energy in the CM. frame, and Eg and p* B are the energy and momentum 
of the reconstructed particles in the CM. frame. 

The signal region for B^ 1 — > J /y^K* 1 candidates is defined as M^ c > 5.27 GeV/c 2 and -35 
MeV < AE < 30 MeV. The final T] J / y invariant mass distribution is shown in Fig. [3] together with 
the fitted results. No hint of a narrow resonance is evident from the current statistics. No X(3872) 
signal is seen and we obtain the limit 3S{B ± -)• X(3872)X ± )^'(X(3872) ->■ J/yq) < 3 - 8 x 10 ~ 6 
at 90% C.L. 




(GeV/cf) 



Figure 3: The fit to the Mjiy^ distribution. Red dashed (green long dashed) curve shows the signal for 
2P — > ^{—^ J I y/T])F^ (phase space component 2P —> J / i//tjA' ± ), while black dashed-dotted curve shows 
the background parameterized using B — > J/\j/X MC sample. 



Complimentary to B — > {j/yi\)K^ study, the search for the X(3872)'s C-odd partner is also 
carried in B — > (XciY)^ process. Besides the X(3872)'s C-odd partner, we also keep an eye on 
any other possible narrow charmonium or charmonium-like candidate. 

After all the event selections, Fig. [| shows the M Xcl y distribution with M^ c > 5 .27 GeV / c 2 (top 
left), enlarged M Xci7 distribution (top right), M bc distribution with 3.66 < M XclY < 3.708 GeV/c 2 
(bottom left), and M bc distribution with 3.805 < M Xcl7 < 3.845 GeV/c 2 (bottom right). The dots 
with error bars are data, and the blue solid lines are the proj ections from 2D unbinned maximum 
likelihood fits. NoX(3872) signal is seen and the limit ^(fl* ^X(3872)X ± )^(X(3872)^y;t: f i)< 
2.0 x 10~ 6 at 90% C.L. is obtained. Besides the clear y/ signal, we find a clear evidence of a nar- 
row peak at 3823 MeV in M Xcl7 , as shown in Fig. ^| with pink dashed curve (named as X(3823)). 
The signal significance is 4.2 a with systematic error included. The mass and width of this peak are 
estimated to be 3823.5 ±2.1 MeV/c" and 4± 6 MeV, respectively. We noticed that charmonium 
model predicts a narrow state ( 3 D 2 cc) at around 3810-3840 MeV/c 2 f|. So theX(3823) is proba- 
bly \jf 2 state. The measured 3S{B ± ^X(3823)X ± )^(X(3823) -> y% c \) is (9.70±|J|±J;g|) x 10~ 6 , 
where the first errors are statistical and the second systematic. 

4. Bottomonium and bottomonium-like states 

The spin-singlet states hb(nP) and rjb(nS) alone provide information concerning the spin-spin 
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Figure 4: The projections from 2D unbinned maximum likelihood fits to the M Xci y distribution with M\, c > 
5.27 GeV/c 2 (top left), enlarged M Xcl y distribution (top right, yellow arrow shows the X(3872) position), 
M bc distribution with 3 .66 < M Xcl y < 3 .708 GeV/c 2 (bottom left), and M bc distribution with 3 .805 < M Xcl y < 
3.845 GeV/c 2 (bottom right). The dots with error bars are data. The blue solid line is form the overall fit. The 
red large-dashed and pink dashed curves are for y/' an d V2 signals, respectively, while black dotted-dashed, 
dark green two dotted-dashed, and cyan three dotted-dashed curves are for the combinatorial background, 
B — > (//'(other than % c \ y)K events, and peaking background component, respectively. 



(or hyperfine) interaction in bottomonium. Measurements of the hf,{nP) masses provide unique 
access to the P-wave hyperfine splitting, AMhf = (M(n 3 Pj)) — M(n l P\), the difference between 
the spin- weighted average mass of the P-wave triplet states {%bj{nP) or n 3 Pj) and that of the cor- 
responding hb(nP), or n l P\. We use a 121.4fb _1 data sample collected near the peak of the T(5S) 
resonance (y / ?~ 10.865GeV) with the Belle detector to report the first observation of the hb(lP) 
and hb(2P) produced via e + e~ — > hb(nP)7i + 7t~ in the T(5S) region [^]. 

We observe the hb(nP) states in the % + %~ missing mass spectrum of hadronic events. The 
missing mass is defined as M^ ss = (.fV(5S) ~ Piz+%-) 2 > where Pj^s) i s ^ e 4-momentum of 
the T(5S) determined from the beam momenta and P K + n - is the 4-momentum of the K + %~ system. 
The n + K~ transitions between T(nS) states provide high-statistics reference signals. 

To reconstruct the T(55) — > hb(nP)n + 7i~ transitions inclusively, we use a general hadronic 
event selection. The M m ; ss spectrum after subtraction of both the combinatoric and Kg — > n + n~ 
contributions is shown with the fitted signal functions overlaid in Fig. [|. The significances of the 
hb(lP) and hb{2P) signals, with systematic uncertainties accounted for, are 5.5a and 11.2(7, re- 
spectively. The measured masses of h b (\P) and h b (2P) areM = (9898.21 J ;^};?)MeV/c 2 andM = 
( 10259.8 ±0.6+};o)MeV/c 2 , respectively. Using the world average masses of the %bj(nP) states, 
we determine the hyperfine splittings to be AM H f = (+1.7±1.5)MeV/c 2 and (+0.5+j'2)MeV/c 2 , 
respectively, where statistical and systematic uncertainties are combined in quadrature. 

The observation of anomalously high rates for T(55) — > Y(nS)n + 7i~ in = 1,2,3) and Y(5>S) — > 
hb(nP)n + n~ (m = 1,2) transitions suggests that exotic mechanisms are contributing to T(55) de- 
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Figure 5: The inclusive M m i ss spectrum with the combinatoric background and contribution subtracted 
(points with errors) and signal component of the fit function overlaid (smooth curve). The vertical lines 
indicate boundaries of the fit regions. 
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cays. Amplitude analyses of the three-body T(55) — > T(nS)K + K decays with T(nS) 
are performed by means of unbinned maximum likelihood fits to two-dimensional M 2 [T(nS)n 



vs. M 2 {Y(nS)n ] Dalitz distributions [10]. One-dimensional invariant mass projections for events 
in the T(nS) signal regions are shown in Fig. ||, where two peaks are observed in the T(nS)n sys- 
tem near 10.61 GeV/c 2 and 10.65GeV/c 2 (named as Z fo (10610) and Z fo (10650)). The combined 
statistical significance of the two peaks exceeds 10a for all Y(nS)7i + 7t~ channels. 
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Figure 6: Comparison of fit results (open histogram) with experimental data (points with error bars) for 
events in the (a) T(15), (c) T(25), and (e) T(35) signal regions. The hatched histogram shows the background 
component. 



To study the resonant substructure of the T(55) — > hb{nP)% + K~~ (m = 1,2) decays we measure 
their yield as a function of the h^lP)^ invariant mass. The decays are reconstructed inclusively 
using the missing mass of the % + %~ pair, Mmiss- We fit the M m i ss spectra in bins of ^(IP)^ 
invariant mass, defined as the missing mass of the opposite sign pion, M m i ss (n T ). We combine the 
Mmiss spectra for the corresponding M m i ss (n + ) and M m i ss (7T~) bins and we use half of the available 
M m i ss (ft) range to avoid double counting. 

The results for the yield of T(55) — > hb(nP)7i + 7i~ (m = 1,2) decays as a function of the 
Mmiss (ft) are shown in Fig. [7|, where the fit results are shown as solid histograms. The two-peak 
structures are clear in both of them. The default fit hypothesis is favored over the phase-space fit 
hypothesis at the 18 <T [6.7(7] level for the h b {\P) [h b (2P)]. 
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Weighted averages over all five channels give M = 1 0607 .2 ±2.0 MeV /c 2 ,T= 18.4±2.4 MeV 
fortheZ fe (10610) and M = 10652.2 ± 1.5MeV/c 2 , T = 11.5 ±2.2 MeV for the Z^(10650), where 
statistical and systematic errors are added in quadrature. Angular analysis favors a J p = 1 + as- 
signment for both Z^ states, which must also have negative G-parity. Transitions through to 
the hb(nP) saturate the observed n + K~hb(nP) cross sections. The two masses of Z£ states are just 
a few MeV above the B*B and B*B* thresholds, respectively. The Z£ cannot be simple mesons 
because they are charged and have bb content. 




Figure 7: The (a) hb(lP) and (b) hf,{2P) yields as a function of M m ; ss (7r) (points with error bars) and results 
of the fit (histogram). 



5. Summary 

We reviewed here some recent results on exotic states at the Belle experiment, including the 
measurement of the cross sections of 77—)- cotp, (pip, and coco; the search for the X(3872)'s C-odd 
partner in riJ/xff and y%c\ modes in B decays; the evidence of 1//2 m ^ ± — > K^YXcll tne observation 
of h b (lP), h b {2P), T(1D), and two charged Z fo (10610) and Z 6 ( 10650) states. 
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